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THE CONVEXITY OF THE FREE BOUNDARY FOR A 
PARABOLIC FREE BOUNDARY PROBLEM 

HSU AN KU LIU 



O 

Cn ' Abstract. In this paper, we study a parabolic free boundary problem which 

5^ ' shows that the solutions of this free boundary problem are increasing functions. 

^ I Furthermore, we provide a rigorous verification for that the free boundary for 

this problem is concave. As an application to the American option pricing 
problem, our results can be extended to consider the convexity of the early 
^^ ' exercise boundary of an American call. 



v«J . 1. Introduction 

,^ ' In this paper, we consider the one phase one-dimensional free boundary problem 

^T^ I for linear parabolic equations. Precisely, given a positive real-valued function ijj e 

qh. C3([0,oo)) with 

(1.1) Coi^i^) I > f- ^ - < '^' 

^ . I < for d < a; < oo, 

^f\ • for some d > 0, we shall discuss some properties of the solution {s,u) for 

m 
in 

^' 

o 
en 



(1.2) 


Cu = 


0<x < s{t), 


< i < oo 


(1.3) 


u{x,t) > -0(x) 


0<x < s{t), 


< i < oo 


(1.4) 


u(0,t) = v(o) 


< i < oo. 




(1.5) 


u{x,0) ^ ^/j{x) 


< a; < s(0). 




(1.6) 


u(s(t),i)=V(s(t)) 


< i < oo. 




(1.7) 


^{s{t),t)^ij\s{t)) 


< i < oo, 




where 







(P) 



and £o is defined as 

We assume that coefficients a,b,c e C^+"([0, 00)) for some a S (0, 1) with a{x) > 
Co > for < X < 00 and c{x) < for < a; < 00. This problem had been studied 
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by Kotlow jT2]. He showed that the problem (P) is well-posed and that u{x,t) and 
s{t) are both nondecreasing functions of t. When C is given as 

C - ^" 

Problem (P) is called a Stefan problem. Friedman [8] showed that the free boundary 
for the Stefan problem is smooth and monotone increasing. Moreover, Friedman 
and Jensen [5] provided conditions under which the free boundary is a concave 
function by using the maximum principle to control the level curve. 

By using the Feynman-Kac formula, the operator C in Problem (P) has a con- 
nection with the option pricing problem when the stock price is measured by time- 
homogeneous diffusions. Black and Scholes (BS) [1] suggested that the stock price 
can be measured by a geometric Brownian motion and derived a second order linear 
parabolic operator, namely BS pricing operator, for pricing an European option. 
In the BS framework, Merton [15] showed that the price of the American option 
together with an early exercise boundary is a solution of Problem (P). 

Over the past two decades, a great deal of effort has been made on solving 
Problem (P) for pricing an American option. However, the exact solution of the 
American pricing model has not been found yet. Several researchers such as [5], [3], 
[6] , [To] , [11] , [13] concentrated on finding more accurately expansions or simulations 
for the early exercise boundary. By observing the results of simulations, people are 
convinced that s{t) is convex for the American put and is concave for the American 
call. The rigorous verifications of the convexity for s{t) of the American put have 
been proposed by Chen et al. [1] and Ekstrom [S] during 2004. The convexity of 
s for the American put provides a useful information for its approximation. When 
the remaining time is close to zero, Chen et al. [4] used this information to provide 
an asymptotic formula for the early exercise boundary of the American put. 

To show the convexity of s{t) for the American put, Chen et al. [4] and Ekstrom 
[5] defined s{t) = Ke^^^^ and showed that z{t) is a convex function by using the 
same argument of [S], where t — T — t. Consequently, s{t) = Ke^^'^^ is a convex 
function since e^ is a convex function. However, their methods do not work for 
verifying the concavity of s(t) for the American call. Although their methods can 
show that z{t) is a concave function for the American call, we still do not know 
the concavity of s{t) = Ke^^'^\ So far, we have not seen any published works on 
the concavity of s{t) for the American call. 

In this paper, we make the following assumptions which will be used later. 
Assumptions 

(A) c{x) is a nonincreasing function in [0,oo). 

(B) il>{x) is a positive strictly increasing convex function in [0, c»). 

(C) £Coi}{x) < and c(x) -I- b'{x) < in [0, oo). 
(D)limsup^^^/:oV'(e)<0. 

Under the given assumptions, we shall show that s(t), u{-,t) and u{x,-) are 
increasing functions . The main contribution of this paper is to provide a rigorous 
verification of the concavity of s(i). The obtained results can be applied to consider 
the convexity of the early exercise boundary in American option pricing problems. 

The content of this paper is organized as follows. In section 2, we provide some 
useful properties for the solution of Problem (P) . The concavity of the free boundary 
is obtained in section 3. Finally, a concise conclusion is given in section 4. 



the convexity of the free boundary for a parabolic free boundary problem 

2. Properties of the Solution 
Let {s,u} be the solution of (P) and denote C, namely continuation region, as 

(2.1) C = {{x, t);0<x < s{t), 0<t<oo}. 
Now we define u on Q by 

u(x t) ^ b'^''--^^ if(a:,t)eC, 

\i^{x) iiix,t)eQ-C, 

where Q — (0, cx)) x (0, oo). 

In this section, we will show that 

(1) s{t) is a strictly increasing function with s(0) = d, 

(2) u{x, •) is a strictly increasing function in C, 

(3) u{-,t) is an increasing function in C. 

In order to prove these results, we need some preliminaries. 

Definition 2.1. {s,u} is a solution of pr^ - pTT)) if 

(i) s e C'''([0,oo)) for some 7 £ [1/2,1) and is Lipschitz continuous for < 

t < CO. 

(ii) ueC^'\C)r\C°iC). 

(iii) ^ is bounded and has a continuous extension to C. 

(iv) 4jT is bounded and has a continuous extension to C — (0, 0) — (s(0), 0). 
(v) Conditions (|1.3p - (|1.7l) are satisfied. 

Definition 2.2. Given t G [0, 00), the t-section of C is defined to be 

(2.2) Ct ^ {x eR\0 <x < s(t)}. 



Clearly, we have 



c = u (a X {t}) 



t<oo 



and 

(2.3) s{t) = sup{x|x e Ct}. 

The following properties of u{x,t) and s{t) have been proved by Kotlow [T^ . 

Lemma 2.3. Let {s, u} be a solution of (P). They have the following properties: 

(a) ut > in C. 

(b) s(0) = d and s{t) > d for < t < 00. 

(c) s(t) is a nondecreasing function. 

(d) There exists a s°° £ (d, 00) such that s{t) — )> s°° uniformly as t ^ 00 if (D) 
holds. 

Lemma 2.4. Let {s,w} he a solution of (P). At points {s(t),t), t > 0, u satisfies 
utis{t),t)^0. 

Proof By (HH]), we have u{s{t),t) = ^is{t)). Differentiating u{s{t),t) = ipislt)) 
with respect to t, we have 

u,{s{t), t)s'{t) + ut{s{t), t) = ^'{s{t))s'{t) 

for each t > 0, where s{t) is differentiable. Since Ux{s{t),t) — 'tp'{s{{t)) by (|1.7p . 
we have ut{s{t),t) = for almost every i > 0. By (iv) of Definition 2.1, we have 

utis{t),t)=0. D 
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This proof is proposed by Kotlow IT^ . 

Lemma 2.5. Let {s,u} be a solution of (P) and w = ii — tp in Q. Then w{-,t) has 
a local maximum in (0, s(t)). Moreover, this local maximum can not lie in {d,s{t)). 

Proof. By lOD and (US]), we have ^(O,^) = w{s{t),t) = and w{x,t) > on C. 
This iniphes that there exists a d' G (O,s(to)) for some to > 0, which may depend 
on to, such that w(d',to) is a local maximum. Now, we claim that d' ^ (d, s(fo))- 
Suppose that d' G (d, s(to)) is a local maximum of w{x,to). We define 

^0 = {{x, t)eC\d<x< s{t), t < to}, 

dpQto ^{ix,t)ednt,\t<to}. 

and apply the differential operator £ to w. By ()l.ip and (|1.2p . w satisfies the 
parabolic differential equation 

a 

£w = {Co - T^)u - Co'ipix) = -Cotljix) > on ^tg- 

Now, we apply the maximum principle to this equation. Since (d' , to) G ^to ~ dpil^to 
is a nonnegative maximum over flta , it implies that w is a constant function on r2t„ . 
This contradicts to that Wt — Ut > 0. So d' < d. D 

Definition 2.6. Let 

(2.4) Xo = a(^)^ + (^(a;) + a'(^))^ + (c(^) + &'(a;)) 

be an elliptic operator. We define a parabolic operator M. as 

M=Mo-?-. 
at 

Theorem 2.7. Let {s,u} he a solution of (P). Then 

(a) s{t) is a strictly increasing function. 

(b) u^ix,t) > for {x,t) eC if (A), (B) hold. 

(c) Ma;(a;,t) < %p'{x) for (x,t) e Cd, w/iere Cd = {(x,t e C|a; > d)}, if (A), 
(B), (C) hold. 

Proof. For (a), we only need to show that s(t2) ^ s(ti), for t2 < ti by (c) in 
Lemma 12.31 Suppose that there is an interval [ti,t2] such that s{t) = s(ti) for 
all t g [ti,t2], then Ux{s{t),t) is a constant function for all t G (ti,t2)- Since 
£u = in (0, s(fi)) X (ti,t2) and u{s{t),t) — i/'(s(^i)) for all t e [ti,t2], we have 
u e C°°([0,s(ti)) X (ti,t2)). Since Ut(s(ti),ti) = for t G (ti,t2), u* > in 
(0,s(ti)) X (ti,t2) and Cut = in (0, s(fi)) x (ti,t2), we have Utx{s{t),t) < for 
t G (ti,t2) by applying the boundary point form of the maximum principle. And 
then, we have Uxt{s{t),t) = Utx{s(t),t) < which implies that Ux{s{t),t) is strictly 
decreasing for t G (ti,t2)- On the other hand, s{t) is a nondecreasing function in 
[ti,t2] and ij{x) is a strictly increasing convex function. This implies that ■0'(s(t)) 
is a nondecreasing function in [ti, t2]. Hence, we obtain that Ux{s{t),t) — ip'{s{t)) is 
a nondecreasing function in [ti,t2] and there is a contradiction. So s{t) is a strictly 
increasing function. 
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We now show that (b) Ux > 0- We first consider that 
d d du d 

= ■7^{a{x)uxx + b{x)ux + c{x)u) 

= a(x)-—^Ux + {h{x) + a'{x))—Ux + (c(x) + h' {x))ux + c {x)u 
Consequently, Ux satisfies the parabohc difi'erential equation 

Mux ~ — c {x)u. 

Let w{x,t) — u{x,t) — ip{x) on Q. Then w{0,t) — and w{x,t) > on C by (|1.3p 
and (|1.4p . So we have 

/„ X ,. w(h,t) — w(0,t) 
Wx(0,t) ^ hni ^ ' , ^^^>0. 

h^0+ h 

This imphes that 

(2.5) M,(0,i)>^'(0)>0 

for < t < CO. We apply the maximum principle to this equation. Since 
Uxisit),t) ^ ^'is{t)) > for i > by ^Ji), UxiO,t) > for < i < oo by 
(|23)) . Uxix,0) = ip'ix) > ioT < X < s(0) by (fTSl) and c'(a;) < by assumption 
(A), this shows Ux{x, i) > for < a; < s{t) and < t < oo. 

Finally, we show that (c) Ux{x,t) < 'ip'{x) ioi d < x < s{t) and < t < oo. Let 
w{x, t) = u{x, t) — tp{x) on Q. Since w{s{t),t) — Wx{s{t),t) = by (|1.6p and (|1.7p . 
w{x,t) > on C by (jl.3p . and the continuity of w and Wx, we have that, for any 
i > 0, there exists a (5 > such that Wx{x, i) < for s{t) — S < x < s{t) and i > 0. 
Suppose that there is a x' £ {d,s{t) — S) with Wx{x',t) > 0, then there exists a 
local maximum in [x',s{t)) C (d, s(t)). This contradicts Lemma 12.51 So we obtain 
that there is no a; G (d, s(t)) with Wx{x,t) > 0. Hence, we have Wx{x,t) < for 
{x,t) G Cd- This implies that Ux{x,t) < 'ip'{x) for {x,t) £ Cd- □ 

3. Concavity of the free boundary 

Suppose that (A), (B) and (C) hold. We will show that s{t) in (j2.ip is a concave 
function; consequently the continuation region of C in (j2.ip is a convex set. 
In order to prove our main theorem, we need the following lemma. 

Lemma 3.1. Let {s,u} be a solution of (P) and define w — u — ip on Cd, where 
Cd = {{x,t) e R'^\d < X < s{t),0 <t< oo}. At points {s{t),t), t > 0, w has the 
following properties 

w{x,t) = 0, wt{s{t),t) = 0, Wx{s{t),t) = 0, 

lim w;^^(a;,i) = ——{Coip{s{t))), 

x-i-s(t) a{s(t)) 

\m\ wtx{x,t) = {Co^{s{t)))s'{t), 

x^s{t) a[s[t)) 

Ihn wttix,t) = --^{Co^{s{t))){s'f{t), 
where s(t) is differentiable. 
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Proof. Since u{s{t),t) — ■i(j{s{t)), Ux{s{t),t) ^ '([;' {s{t)) and ut{s{t),t) = 0, we obtain 
the first three equalities. By (jl.ip and ()1.2p . we have 

o-Wxx + bwx + cw — Wt ~ —Coip. 

Since w = 0, Wx = and wt = at (s(i), t), we have 

Wxx = --{Cotp). 
a 



wx{s{t),t)^0 



Differentiating the equality 
with respect to t, we have 

WxxS'{t) +Wtx = 0. 

So Wtx = —Wxxs'{t) — ^{Co'4>{x))s' {t) . Next, Wu can be found by differentiating 

wt{s{t),t)=0. 

We obtain 

WxMt),t)s'{t)+wu{s{t),t) = Q. 
So Wu = ~wxt{s{t),t)s'{t) = -\{Coij{x)){s'f{t). U 

Now, we can state and proof our main theorem of this paper. 

Theorem 3.2. Let (A), (B), (C) and (D) hold, Mgip' < and {s,u} be a solution 
of (P). Then s{t) is a concave function. 

Proof. We have known that s{t) is a strictly increasing function. Suppose that 
there is an interval / such that s{t) is a convex function in the interval /. There 
is a io G -^ with s'(to) = m > since s(t) is a strictly increasing function and is 
differentiable almost everywhere. Then, s'{t) > m for almost every t > t^ in I. We 
consider the line 

y{t) ^ m{t - to) + s{to) 
for some i > 0. So y{to) = s{to). Since s{t) is bounded above by (d) of Lemma 
12.31 and ?7i > 0, there must exist another point ii > to such that j/(ii) = s{ti). 
Let f{t) = w{y{t),t) for some t > t2 where ^2 — '^'ni{t\{y{t),t) e Cd}- Since 
wt{s{t),t) = and Wx{s{t),t) = 0, we have 

/'(tj) = mwx{y{ti),ti) + wt{y{ti), ti) 

= mwx{s{ti), ti) + wt(s(ij), tj) =0, i = 0, 1. 

We also have /(io) =w{y{to),to) = 0, /(ti) =w{y{ti),ti) = and {y{t),t) e Cd for 
t G {to,ti), which implies that f{t) = w{y{t),t) > for t e (io,ii). So there exists 
a local maximum of / between to and ii, namely /(is) where t^ G {to,ti). This 
implies that /(is) > and /'(^a) — 0. Since w = u — ip is a solution of parabolic 
equation and f{t) = w{y{t),t), which does not oscillate as i — > ti. This implies 
that f{t) do not produce an infinite sequence of local maximum whose locations 
tends to ti. Consequently there is no infinite sequence of local maximum whose 
locations tends to ii. Therefore, we can assume that ^3 is the first one from ii and 
there is no local maximum between ^3 and ii. Since f{to) = f{ti) — 0, /(i3) > 0, 
and f'{ti) = 0, i = 0, 1, 3, we have 

(3.1) /'(t)<0fortG(i3,ii) 

and f'{t) > for some interval, say (^4,^3), where to < ^4 < ^3- 
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Since 

lim ^^ = -.'(i) 

by THopital's rule and w > on Cd, let 

f^ iiix,t)eCd, 

\-s'{t) iix = s{t), 

which is well-defineci on Cd- Then we have that 

..^^ fit) ^mw^{y{t),t)+wt{y{t),t) 

^ ' ^w,{y{t),t){m + v{y{t),t)) 

for t > t2- Applying the differential operator C to the equality vwx = Wt, we obtain 
that V satisfies the following parabolic equation 

(3.3) avxx + (foH —)vx {Mo4'')v - vt = on Cd- 

Wx Wx 

Let Ta be the level curves on which v = a- Since w^ < 0, M.o^' < and v 
satisfies the parabolic equation p.3p . the x-coordinate along Fq, can not first in- 
crease (decrease) and then decrease (increase) by using the extensions of maximum 
principle (see Appendix A). This is because that there would be a region whose 
parabolic boundary is a part of F^; consequently v = a in this region and hence 
V = a in Cd- Since F^ is understood to be continued as long it remains in Cd, for 
each a there is a gait) such that 

Ta = {igait),t)\vigait),t) =a,t> 0}. 

Since f'iti) — 0, i — 0,1,3 and f'{t) = Wxiyit),t){m + v{y{t),t)), we have 
viyiti),ti) — — m, i = 0, 1,3. This implies that {y{ti),ti) G T^„i, i = 0, 1,3. Now 
we consider the function g_„i(i). Since (y(ti),ti) G F_„j, that is viy{ti),ti) = — m, 
i = 0, 1, 3, we have y{U) = g^^iti), i = 0,1, 3. 

Since g-mit) is continuous on (t2,ii), we only have the following two cases: (1) 
yit) < g-mit) for t G (^3,^1), and (2) y{t) > 5_,„(i) for t G (i3,ii)- 

We first consider case (1). We have Wx < on Cd by (c) of Theorem 12.71 Since 
f'{t) = Wxiyit),t){m + viyit),t)) < for t G (i3,ti) by dXT]) and (02]), this implies 
that 

(3.4) viyit),t) > ~m, for i G (^3,^1)- 

Since g^miti) = viti) = s(ti) and yit) < g-mit) < s{t) for t G (^3,^1), there is a 
5 > such that y'{t) > g'{t) > s'{t) for t e {h - 6,ti)- Since y'{t) = m, we have 
v{s{t),t) = -s'(t) > -y'it) = -TO for t G (ti - 5,ti)- Let fl = {{x,t)\y{t) < x < 
sit), ti — S < t < ti}- Since v{g-miT),T) = —to, for r G (t, ri + S), the negative 
minimum of w in 17 which is attained at [x' ,t') is less than —to. By Theorem lA.il 
which is displayed in Appendix A, we have 

v{x,t') = v{x' ,t') < —TO, Va; such that ix,t) G fi, 

which contradicts to u(s(i'),t') > —to and v{yit'),t') > —to. 

Now we consider case (2). We have known that the level curves Fq of a parabolic 
equation is continuous. We consider the line y{t) for t G (^3,^1) U (^2,^0)- Since 
p.4p is also true for case (2), we have v{y{t),t) > —to for t G (^3, ti). And then, we 
also have /(io) = and /(t) > for t G (t2, ^o)- This implies that there is a <^2 > 
such that f'{t) < for t G ito~S2,to)- Since fit) = Wx{y{t) ,t){m + v{y{t) ,t)) and 
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f'{t) <Oforie (to— S2, to), we have v{y{t),t) > —miorte {to— 62, to). Now we can 
select a suitable (5 > such that w (2/ (t),i) > — to for f e {to — S,to)D{ti—S,ti). Since 
v{y{to),to) = -m = v{y{ti),ti) and v{y{t),t) > -to for t e (io - (5, ^o) U (ti -(5,ii), 
there exists a t' G (to ~ <^i ^o) and a i" e (ti — (5, ti) such that 

v[y[t'),t') = /? = v[y[t"),t"), for aomc ^ > -to. 

Since the level curves of a parabolic equation is continuous, there is a level curve T p 
connecting {y{t'),t') and ly(t"),t"). On the other hand, we have {y{t"),t") G fii, 
where Q.i = {{x,t) G C\g-m{t) <x< s{t), to < t < ti}. This contradicts to that 
r_,„ n r^ 7^ 0. So case (2) does not hold. 

Since both case (1) and case (2) do not hold, we conclude that there is no such 
point to with s"{t) > for f > 0. Thus, s{t) is a concave function. D 

Remark 3.3. Given a Cz R and gait) be the function such that 

v{ga{t),t) = a 
with ga(to) = s{to)^ where ti(s(to),io) — ct- Then 

dv dga{t) 

By Sard's lemma, the set of Vx{x,t) = has measure zero. Hence —^ is defined 
for almost every point on il. We consider the following IVP 

/o ^\ dga{t) vt 

(3.5) — -— = (a.e.) 

dt Vx 

with ga{to) = s{to)- Indeed, the weak solution for (13. 5p exists. Therefore ga{t) is 
continuous for all t with v{ga{t),t) = a. 

The detail about the extensions of maximum principle can be found in Friendman 

Remark 3.4. Let ip{x) denote the final payoff function and suppose that the stock 
price is measured by the time-homogeneous diffusions 

dX{t) = r{X{t))dt + a{X{t))dB{t), and X{s) = x. 

By using the Feynman-Kac formula, we have that the European option's price v 
under the constant interest rate r > is a solution of the following initial value 
problem 

Cv = 0, with v{x,0) = ip{x), 

where t = T — t, £ = £0 — -^ and 

(3.6) ^"-r'(-)^+'^(-)|:-^- 

However, an American option gives the holder a right to exercise prior the ma- 
turity. For the possibility of early exercise, the value of an American option is no 
longer less then its immediate exercise value, that is v{x, t) > ip{x). Now, we sepa- 
rate the domain {(x,i)|0 < x < 00, < t < T} into two parts: (i) a continuation 
region C = {{x,t)\v{x,t) > ip{x)}, and (ii) a stopping region S = {{x,t)\v{x,t) = 
'ip{x)}. Given t G (0, T), when '4'{x) is a strictly increasing function, there is a time- 
depended function s{t) such that v{x,t) > ipix) for < a; < s{t) (see [3]). Hence, 
the continuation region is described as C = {(a:,t)|0 < a; < s{t), < t < T}. On 
the other hand, given t G (0,T), when ipix) is a strictly decreasing function, there 
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is a time-depended function s{t) such that v{x, i) > 'tp{x) for s{t) < x < oo. Hence, 
the continuation region is described as C = {{x,t)\s{t) < x < oo, < i < T}. And 
then, by using the no arbitrage condition, the price of the American option satis- 
fies the high contact condition, that is Vx{x^t) = ip'ix). Therefore, the American 
option's price v with an early exercise boundary s is a solution of Problem (P), 
where the operator Cq is given by p.6|) . Consequently, the rigorous verification 
of the concavity for the early exercise boundary can be obtained by applying the 
presented results. 

4. Conclusion 

In this paper, we consider a parabolic free boundary problem. Under given as- 
sumptions, we have shown that the solutions of this problem are strictly increasing. 
Moreover, we have obtained that the free boundary s{t) is a concave function. The 
American call pricing model derived by Merton [15] is a specific case of this para- 
bolic free boundary problem. Nevertheless, the rigorous verification of the concavity 
for the early exercise boundary has been neglected by their works. By using our 
results, a rigorous verification of the concavity of the early exercise boundary can 
be obtained. 

Appendix A. Extensions of maximum principle 

In this section, we display the extensions of maximum principle in [7] . Let D he a 
2-dimensional domain bounded by an open interval i? on t = 0, an open interval Bt 
on i = T and two continuous curves Ci : x = ji{t), i — 1^2 defined for < i < T. 
For any P° — {x^,t°) in D, we denote C(P°) as the component oi D (1 {t ^ t^} 
which contains P°. 

Consider the differential operator 

d'^ d d 

with continuous coefficients in domain D. Assume that (A) C is parabolic in D; 
(B) the coefficients of L are continuous functions in D\ (C) C{x) < in _D. Here, 
Friedman [7] did not make the assumption of d{x). 

Theorem A.l. Let (A), (B) and (C) hold. If Cu > (Cu > 0) in D and if 
u has a positive maximum (negative minimum) in D which is attained at a point 
po = (x°,i°), then u{P) = u{P°) for all P G C{P°). 
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